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SOME  EXPERIMENTAL  ARTIFACTS  IN  LEED 
STUDIES  OF  PHASE  TRANSITIONS 

ABSTRACT 


Low  energy  electron  diffraction,  LEED,  is  an  important  experimental  technique 
for  the  study  of  continuous  order-disorder  phase  transitions.  Extraction  of  the 
critical  exponents  and  their  associated  amplitudes  requires  an  accurate 
determination  of  the  single-scattering  intensity  vs.  parallel  component  of  the 
scattering  vector.  Interpretation  of  measured  intensities  is  complicated  by 
multiple  scattering,  variation  of  the  scattering  factor,  and  thermal  vibrations. 
These  can  cause  the  measured  intensity  profiles  to  be  skewed  and  can  add 
substantial  intensity  which  is  peaked  at  the  Briilouin  zone  center.  An  example  of 
subsidiary  structure  induced  by  multiple  scattering  is  given  for  Au(llO).  The  effects 
on  the  measured  critical  exponents  13,  y  am*  v  are  shown  to  be  small  in  most  cases. 
However,  the  amplitudes  associated  with  these  exponents,  such  as  the  correlation 
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Surfaces  provide  models  of  two-dimensional  systems  which  have  a  rich  variety  of 
continuous,  second-order  phase  transitions.  There  has  been  a  recent,  rapid  increase 
in  the  use  of  low-energy  electron  diffraction,  LEED,  to  study  these  transitions  and 
the  results  of  these  studies  are  beginning  to  provide  tests  of  the  concepts  of 
Universality  and  Scaling  (lj.  For  example,  the  Au(llO)  (1x2)  order-disorder  transition 
is  now  known  to  belong  to  the  Ising  Universality  Class  (2,3).  A  number  of  factors 
which  influence  the  LEED  intensities  have  not  been  accounted  for  in  the  analyses 
reported  so  far.  This  paper  discusses  two  of  these  factors:  multiple  scattering  and 
thermal  diffuse  scattering,  TDS.  These  effects  are  not  new  (4,5)  and  the  primary 
purpose  of  this  paper  is  to  discuss  how  they  influence  LEED  studies  of  critical 
behavior. 

The  data  required  for  analysis  of  critical  properties  consists  of  the  intensity 
measured  along  a  direction  parallel  to  the  surface.  The  variation  of  such  a  'spot 
profile"  for  temperatures  near  the  critical  temperataure  Tc  contains  the 

information  needed  to  characterize  the  phase  transition.  For  an  continuous 
transition,  the  intensity  in  the  profile  is  usually  written  in  the  form  (1) 

l(Sj|,T)  -  l0  (T)  «  (S|j  -  Ghk)  +  X<S||.  T)  (I) 


where  lQ(T)  is  the  intensity  due  to  long-range-order  in  the  ordered  phase  and 
X(5||,T)  is  due  to  fluctuations  occuring  near  Tc;  Ghk  locates  the  reciprocal  lattice 
rods,  5||  is  the  component  of  the  scattering  vector  parallel  to  the  surface;  5|j  = 
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k|j-k0||  where  ky  and  kQj|  are  the  surface  components  of  the  scattered  and  Incident 

wave  vectors  respectively.  Near  Tc,  the  quantities  lQ,  x  and  t  (the  Inverse  of  the 
width  of  x)  are  predicted  to  have  the  following  power  law  behaviorf  1-3,5): 


10  —  Al  1 12B 

(2a) 

x(o,t)  ^  x0itry 

(2b) 

4(t)  -  ytf° 

(2c) 

where  13,  y  and  u  are  the’ critical  exponents  whose  values  determine  the  universality 
class  of  the  transition.  £  is  also  known  as  the  correlation  length. 

Diffraction  from  a  real  surface  is  more  complex  than  indicated  by  Equation  1, 
which  neglects  the  angular  and  energy  dependence  of  the  differential  scattering  " 
cross-section,  the  effects  of  multiple  diffraction,  and  the  effects  of  thermal 
vibrations.  Each  of  these  quantities  depends  explicitly  on  the  scattering  geometry. 

The  latter  two  are  also  functions  of  temperature. 

Finite  size  effects  (6)  due  to  extended  or  point  defects,  such  as  steps  and 
impurities,  can  also  influence  the  diffracted  intensity  but  are  not  considered  here. 

The  important  effects  of  instrumental  parameters,  some  of  which  can  be 
conveniently  expressed  in  terms  of  a  transfer  function  (7),  are  also  not  discussed 
here. 

1.  Experimental  Geometries 

Measurements  of  intensity  profiles  are  usually  carried  out  in  one  of  the  three 

Fic^aV 

geometrical  configurations  shown  in  Figure  1.  The  constant  energy  configuration,  ^ 
Figure  la,  fixes  kQ  and  the  direction  of  k  is  varied  at  constant  E.  This  configuration 
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describes  both  the  diffraction  pattern  observed  visually  on  the  screen  of  a  retarding 
field  analyser  and  the  intensity  measured  with  a  movable  Faraday  collector.  The 

constant  ^configuration,  Figure  Kb),  fixes  the  angle  between  k  and  k0  and  E  is 

varied.  The  constant  0  configuration  is  appropriate  when  a  spot  photometer  is  used 
to  measure  the  intensity  at  a  fixed  point  on  the  screen  as  a  function  of  energy.  The 
constant  5  configuration,  Figure  1c,  fixes  the  scattering  angle  between  k  and  k0 

and  the  crystal  is  rotated  about  the  incidence  direction.  This  configuration  applies 
when  the  crystal  is  rotated  with  fixed  detector  position;  the  measured  profile  is  a 
type  of  rocking  curve.  In  each  of  the  above  cases,  the  intensity  is  not  measured  as  a 
function  of  Sn  but  along  some  other  trajectory  in  reciprocal  space.  The  magnitude 

of  S|j  can  of  course  be  calculated  from 

S||  =2TT(E/150.4)1'2(sintf -sintfQ) 

where  3  is  the  exit  angle  and  3Q  the  incidence  angle,  both  measured  with  respect  to 

the  surface  normal.  The  only  technique  which  measures  I  vs  Sjj  directly  Is 

Mirror-Electron-Microscopy-Low-Energy-Electron  Diffraction,  MEMLEED  (2). 

If  the  intensity  along  the  reciprocal  lattice  rod  has  structure  such  as  that  caused 
by  multiple  scattering,  each  of  the  experimental  configurations  in  Figure  1  will  lead 
to  an  incorrect  determination  of  the  angular  profile.  This  is  illustrated 
schematically  in  Figure  2.  The  upper  panel  shows  intensity  contours  around  an 

T  uo  o 

intensity  maximum  on  a  reciprocal  lattice  rod.  The  dashed  line  represents  the 
trajectory  of  the  measurement  and  the  horizontal  line  is  the  desired  trajectory  for 
analysis  using  Equation  I.  In  the  lower  panel,  the  measured  angular  profile  (dotted 
curve)  and  the  true  profile  are  shown.  For  the  case  shown  in  Figure  2,  the  widths  of 
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the  profiles  differ  by  a  factor  of  two.  Near  a  continuous  phase  transition,  the  width 
of  the  line  profile,  Equation  2c,  changes  as  the  temperature  is  varied  near  the 
critical  temperature  Tc.  At  Tc,  the  profile  is  very  narrow  and  the  contours  in  Figure 
2  become  more  elongated  causing  the  error  in  profile  measurement  to  be  reduced. 
Away  from  Tc,  the  contours  become  more  circular  leading  to  a  larger  error  in  profile 

width  measurement.  The  net  result  is  a  temperature  dependent  overestimate  of  4- 
Fortunately,  it  is  simple  to  correct  for  this  effect:  Each  intensity  point  l(S||,  Sj) 

on  the  measurement  trajectory  must  be  normalized  to  the  intensity  measured  on  the 
reciprocal  lattice  rod  at  the  same  Sj_.  To  do  this,  the  inner  potential  must  be  known, 
but  a  several  eV  uncertainty  has  only  a  small  eifect. 

The  discussion  in  Section  2  provides  a  specific  example  for  the  case  of  a  line 
profile  measurement  on  the  Au(l  10)0x2)  surface.  The  uncorreccted  profile  with 
maximum  at  62  eV,  Figure  2a,  is  ten  percent  narrower  than  the  corrected  profile, 
Figure  2c. 

2.  Multiple  Scattering 

It  is  well  established  that  multiple  scattering  is  very  important  in  LEED  [5]  and 
results  in  significant  modulation  of  the  diffracted  intensities.  Most  experimental 
studies  of  phase  transitions  have  attempteJ  to  minimize  multiple  scattering 
within  surface  layer  by  using  normal  inciaence  s;nce  the  differential  scattering 
cross  section  has  its  minimum  value  when  the  scattering  angle  is  near  ninety 
degrees.  However,  multiple  scattering  with  the  substrate  always  occurs  even  in  the 
simplest  two  beam  approximation  (5).  In  the  absence  of  multiple  scattering,  the 
intensity  scattered  from  a  2D  system  would  change  smoothly  with  Sj_,  with  the  only 
variations  due  to  the  differential  scattering  cross  section  and  the  thermal  diffuse 


scattering.  A  glance  at  actual  Intensity  profiles  such  as  that  for  Au(UO)  in  Fig.  3b 
shows  that  significant  modulation  is  the  rule.  Thus,  LEED  can  only  be  a  useful  tool 
for  quantitative  studies  of  phase  transitions  if  multiple  scattering  has  unimportant 
effects  on  the  spot  profiles.  Most  LEED  studies  have  tacitly  assumed  this  to  be  the 
case.  However,  whether  this  is  ever  true  remains  to  be  demonstrated.  In  fact, 
calculations  which  simulate  the  effects  of  disorder  on  the  multiple  scattering  show 
that  the  intensity  profiles  can  be  altered  [8], 

An  example  on  multiple  scattering  effects  is  shown  in  Figure  3  for  the  (0,  -  0.5) 
beam  of  Au(llO)  (1x2).  The  line  profile  shown  in  Figure  3a  was  obtained  in  the 
constant  0  configuration,  Figure  lb.  Instead  of  the  profile  being  a  simple,  nearly 
symmetrical  peak  with  maximum  at  62  eV,  additional  structure  is  apparent  near  45 
eV.  The  center  panel  of  Figure  3  shows  the  intensity  measured  along  the  reciprocal 
lattice  rod.  None  of  the  fine  structure  between  40  and  85  eV  can  be  explained  by  a 
single  scattering  analysis  based  on  the  known  structure  of  this  surface  (9);  it  is  due 
to  multiple  scattering.  The  additional  structure  in  the  spot  profile  is  caused  by  the 
large  peak  at  45  eV  in  Figure  3b.  If  the  assumption  is  made  that  multiple 
scattering  does  not  Influence  line  shapes  measured  as  a  function  of  5||,  then  this 

multiple  scattering  induced  structure  will  be  removed  by  dividing  the  profile 
measurement  of  Figure  3a  by  the  intensity  of  Figure  3b  as  described  in  conjunction 
with  Figure  2.  The  result  of  this  normalization  Is  shown  in  Figure  3c  and  eliminates 
most  but  not  all  of  the  subsidiary  features.  The  remaining  sharp  structure  at  51  eV 
must  be  due  to  breakdown  of  the  assumption  that  the  multiple  scattering  does  not 
vary  with  Sjj.  It  cannot  be  due  to  Incorrectly  accounting  for  thermal  diffuse 

scattering  as  discussed  in  the  next  section  since  failure  to  correct  for  TDS  would 
produce  a  dip  rather  than  a  peak. 

Several  approaches  are  needed  befor  it  will  be  possible  to  better  understand  the 
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relative  importance  of  multiple  scattering  In  Intensity  profile  measurements. 

First,  multiple  scattering  is  a  strong  function  of  geometry.  Thus,  if  spot  profiles 
are  measured  for  a  wide  range  of  experimental  parameters,  the  critical  exponents 
extracted  from  them  will  not  be  the  same  if  multiple  scattering  contributions  are 
significant.  Such  measurements  have  net  been  reported.  The  second  approach  is  to 
compare  LEED  profiles  with  profiles  obtained  by  glancing  incidence  x-ray 
diffraction,  for  which  multiple  scattering  effects  are  not  important.  Finally, 
theoretical  calculations  of  the  effects  of  multiple  scattering  on  spot  profiles  need 
to  be  carried  out  for  both  ordered  and  disordered  systems.  Computational  methods 
needed  for  such  a  study  are  currently  under  development  (8). 

It  should  also  be  pointed  out  that  the  method  of  Bartelt,  Einstein  and  Roelofs 
(3, 10]  for  determination  of  tne  specific  heat  critical  exponent  at  does  not  seem  to 
be  sensitive  to  multple  scattering  effects. 

3  Thermal  Diffuse  Scattering 

The  effects  of  thermal  diffuse  scattering  TD5  due  to  vibrations  of  the  scatterers 
about  their  rigid  lattice  sites  has  also  been  left  out  of  Equation  1.  In  general  TD5 
redistributes  intensity  throughout  the  Brillouin  zone.  Because  the  diffracted 
intensity  from  a  surface  is  distributed  along  a  reciprocal  lattice  rod  rather  than 
isolated  at  reciprocal  lattice  points,  TDS  is  substantially  more  important  for 
diffraction  from  surfaces  than  for  diffraction  from  the  bulk  (1 1). 

For  a  surface  which  does  not  undergo  a  phase  transition,  the  TDS  is  conveniently 
described  in  terms  of  three  components,  each  with  a  different  S  and  T  dependence 
H II).  In  the  high  temperature  limit,  the  zero-phonon,  or  Bragg  intensity  is  the 
rigid  lattice  intensity  reduced  by  the  Debye-Waller  factor: 


where  2MT  =  <(5  •  u)2>  for  single  scattering,  Here  u  is  the  instantaneous 
displacement  of  an  atom  from  its  lattice  site  and  <  . .  >  represents  a  thermal 
average.  The  one-phonon  intensity  is  due  to  scattering  events  which  involve  a 
single  phonon  and  is  peaked  near  the  center  of  the  Brillouin  zone.  Its  T-dependence 
is  given  by: 

l(,)  oc  2MTe_2r1T.  (4) 

Higher  order  phonon  scattering  gives  rise  to  the  multiphonon  intensity  1^  which  is 
nearly  uniformly  distributed  throughout  the  Brillouin  zone  and  has  T-dependence 

|(m)  ex  |_e~2MT_2^Te"2MT>  (5) 

Combination  of  the  analysis  leading  to  Equations  3-5  with  that  leading  to 
Equation  1  plus  the  assumption  that  the  time  scale  for  thermal  vibrations  to  reach 
equilibrium  is  much  shorter  than  the  time  scale  of  the  critical  fluctuations  shows 
that  the  TD5  has  the  same  general  effect  on  the  critical  scattering  x(5 ,T)  as  it  does 
for  the  Bragg  scattering  (3,12).  The  integrations  over  the  phonon  wave  vectors 
required  to  evaluate  these  effects  quantitatively  are  complex  and  cannot  be  carried 
out  analytically.  Furthermore,  the  variation  of  u  near  a  phase  transition  is  not 
known.  This  makes  it  impossible  at  present  to  directly  calculate  the  effect  of  one 
phonon  scattering  on  the  angular  profile  of  x(5,T).  However  it  is  possible  to  draw 
some  qualitative  conclusions. 

First,  the  T-dependence  of  the  TD5  near  Tc  is  slow  compared  to  that  due  to  the 
fluctuations.  Expansion  of  Equations  3-5  in  powers  of  reduced  temperature  t  about 
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Tc  shows  that  the  critical  exponents  will  be  In  error  by  only  a  Tew  parts  per 

thousand  for  realistic  values  of  2M  and  t  <  0.1.  This  is  well  within  the  uncertainty 
with  which  the  exponents  can  currently  be  determined.  However,  such  errors  could 
begin  to  be  important  in  studies  of  both  finite  size  effects  and  higher  order 
correction  terms  to  the  exponents. 

Second,  the  TDS  intensity  always  adds  a  term  to  the  measured  intensity  which  is 
peaked  near  the  Brillouin  zone  center.  Thus,  the  amplitudes  associated  with  each  of 
the  critical  exponents  in  Equation  2  will  be  increased  by  the  addition  of  a 
temperature  dependent  term  if  no  correction  is  made.  For  example,  if  the  TDS 
behavior  measured  for  Ag(lll)  ( 1 2]  is  assumed  be  approximately  correct  for  Au(llO) 
(1x2)  near  Tr  these  effects  are  found  to  be  as  large  as  twenty  percent.  Since  the 

amplitudes  are  not  universal  quantities  but  depend  upon  the  specific  physical  model, 
the  ability  to  experimentally  test  each  model  will  be  limited  unless  the  effects  of 
TDS  are  better  understood. 

4.  Summary 

In  summary,  some  corrections  to  LEED  lineshapp  analysis  have  been  discussed 
which  limit  the  accuracy  witn  which  the  critical  behavior  of  a  continuous, 
order-disorder  phase  transition  can  be  studied.  Multiple  scattering  is  shown  to 
influence  the  measured  lineshape.  Thermal  diffuse  scattering  has  only  a  small 
effect  on  the  extraction  of  critical  exponents  but  causes  significant  errors  in  the 
associated  amplitudes. 

These  results  are  also  expected  to  hold  for  reflection  high  energy  electron 
diffraction,  RHEED,  and  the  results  concerning  thermal  diffuse  scattering  should  be 
valid  for  glancing  incidence  x-ray  diffraction. 
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Figure  1:  Experimental  geometries  for  line  profile  measurements. 

a.  Constant  E:  The  incident  wave  vector  kQ  is  fixed  and  the  scattered  wave 
vector  k  is  moved  along  the  Ewald  sphere  as  shuwn  by  the  dashed  line. 

b.  Constant  0:  The  directions  of  k0  and  k  are  fixed  and  E  is  varied.  The  trip  of 

the  wave  vector  k  moves  along  the  straight  dashes  line. 

c.  Constant  5:  k0  and  k  are  varied  to  move  the  tip  of  S  along  the  dashed  circle 
centered  on  the  origin. 

Figure  2:  The  effect  of  intensity  variations  along  a  reciprocal  lattice  rod  on 
measured  line  profiles.  The  upper  panel  shows  intensity  contours  near  a  maximum 
in  the  intensity.  The  dashed  line  shows  the  trajectory  of  the  measurement;  see 
Figure  1.  The  lower  panel  shows  the  line  profile  measured  by  experiment  (dashed 
line)  and  the  true  line  profile  (solid  line)  measured  along  S||. 

FioureJT  Intensity  measurements  for  Au(!  10)0x2). 

a.  Line  profile  measured  in  the  constant  9  configuration  shown  in  Figure  lb. 

b.  Intensity  vs  energy  profile  measured  along  the  (0.-1/2)  reciprocal  lattice  rod. 

c.  Line  profile  from  a  divided  by  the  intensity  profile  b  at  equal  values  of  Sg_. 


